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Abstract 


We calculated the force to which Cs atoms are subjected in the one-dimensional magneto-optical trap 
(ID-MOT) and properties of this force are also discussed. Several methods to increase the number of Cs 
atoms in the ID-MOT are presented on the basis of the analysis of the capture and escape of Cs atoms in 
the ID-MOT. 

PACs: 3280P; 3270N; 4250 


1. INTRODUCTION 

Laser cooling and trapping of neutral atoms is one of the most active research fields in 
physics in this decade Several kinds of neutral atom traps have been achieved 
experimentally [l * ,one of which is the magneto-optical trap(MOT) whose basic principles 
were discussed in detail in Refs. [2] and [3], The MOT is formed in the intersection of three 
orthogonal pairs of counter propagating laser beams with opposite circular polarization A 
pair of anti-Helmholtz coils is used to generate an inhomogeneous magnetic field whose 
zero point coincides with the center of the intersection region where a stable potential well 
is formed. Atoms confined in this volume will experience a damping force and a restoring 
force. MOT has become almost a standard technique for obtaining large number of cold 
atoms due to its large depth, large trap size and long trap lifetime. A lot of theoretical and 
experimental work has been carried out since Raab et al. 121 first achieved MOT in 1987 In 
1990, Wieman et al.* 4 * built a MOT into which for the first time Cs atoms were loaded 
directly from low pressure Cs vapor in a quartz cell and this greatly promoted the research 
on MOT At present about 10 10 atoms can be trapped in MOT with a density of 10 n 
atoms/cm' 151 and temperatures below the Doppler limit 161 . In recent years, the research of 
MOT has been concentrated on two aspects: ( 1 investigating the properties of MOT with 
the aim of increasing the number and density of trapped atoms and lowering the 
temperature so as to optimize the performance of MOT as a source of cold atoms; (2)using 
MOT to carry out some fundamental or applied research work such as atomic fountain 171 , 
cold atoms collision 181 , atomic interferometry 19 * and Bose-Einstein condensation, etc. 

A simple case of a fictitious atom with a J g = 0 — > J t - 1 transition in the ID-MOT 
*project supported by National Natural Science Foundation of China. 
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In this paper we use a simple method to deal with this problem and have got some 
valuable results. 

In section 2 we calculate the force as well as the capture and escape of Cs atoms in 
ID-MOT ; in section 3 we continue by analyzing the features of the force and we present 
several measures to increase the number of Cs atoms in MOT; finally in section 4 we give 
a conclusion and some comments on our future work. 


2. CALCULATION OF THE FORCE AND THE MOTION OF Cs ATOMS 

A. Calculation of the force 

The force that an atom experiences in a laser field is 110 ’ 


_< ir> 


( 1 ) 


where < > represents the mean value, V = D ■ E is given by the interaction Hamiltonian 
between the induced electric dipole D of the atom $nd the electric field vector E of the 
laser. The authors of Ref. [10] analyzed the force of a fictitious atom with 
J g = 1 -> J e = 2 transition under o + - a~ configuration. The transition used for 
cooling and trapping of Cs atoms is 6 S m ,F g = 4 -> 6P 3a ,F e =5 and the schematic 
diagram of the corresponding energy levels is shown in Fig. I 1 " 1 . 
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Fig. 1 Energy levels of Cs 6 S ll2 >F S = 4 — > 6P 3/2 , i\ = 5 
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There are several kinds of mechanisms which contribute to the total force exerted on Cs 
atoms in the case of complicated energy level structure The most important two are , 
namely , the scattering force and the polarization gradient force 1 ' 1 By extending the 
result of Ref. [10] to the case F g -> F e = F g + 1 ( the contribution of the excited-state 

population is taken into account in the calculation ) , we have 


f f scatter f , i 


grad 


( 2 ) 
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where p { m s) =< F g ,m\p\Fm >,p ( m e) =< F e ,m e \ p\F e ,m e > are diagonal elements 


of the density matrix, S m 


n 2 / 2 


m ± 1 


:1 + r 2 / 4 


is the saturation factor, Q = 2 dEh 


- 1 


is the Rabbi frequency, ( C ”‘) 2 » s corresponding normalized transition probabilities 

shown in Fig.2, mjn±l = 5 + kV + [mg g - (m±l)g e ]p B X— is the relative frequency 

detaining, is the frequency detaining, k is the optical wave vector, v is the atomic 
velocity, p B is the Bohr magneton, g g and g e are the g factors of the ground state and the 
excited state, respectively. , and p mm , 2 = {F g ,m\p\F g ,m - 2)zx$(-2ikVt) represents the 

coherence between sublevels of the ground state. 


In formula (3) , L £ ( p <*> _ )S m m ^ (C;* 1 ) 2 represents the scattering rate for 

2 m = - F g 

r f, 

a + photons by spontaneous emission and _ ^(p ( w f ) - (C" _1 ) 2 is ^at for a 

2 m=-F t 

photons. The atom will obtain a momentum hk or 4ik if it scatters a cf photon or a a 
photon , respectively The net effect is that the atom will experience a scattering force 

f 

I scattcr 


From formula (4) , we know that f^.a is due to the off-diagonal element p m , m . 2 , that is , 
the coherence between sublevels of the ground state. Fig.2 shows the stimulated 
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absorption-stimulated emission process which contributes to the coherence under a* -o ' 
configuration. 



Fig. 2 Stimulated absorption-stimulated emission 
process between sublevels of the ground 
state of Cs under a + -cj ' configuration . 


The authors of Ref [10] also pointed out that f^ would contribute dominantly to the 
force experienced by an atom only when the condition (<? ±1 - 8) < P ( P refers to 

the optical pumping rate. The condition can also be written as AV=P/k , AX=/tP(ge p B 
dB/dx)' 1 ) is satisfied. We may neglect the effect of f p ,d , if this condition is not fulfilled 
Taking the Cs atom for example, with P<0.5MHz under the condition 6= -10 MHz , I = 
4mW/cm 2 , dB/dx = lmT/cm , g g (6S1/2, F = 4) = 0.25 , g. (6P 3/2 , F’ = 5 ) = 0.4 , the 
velocity range is AV<0.45m/s and the spatial range is AX<lmm. Such a small interaction 
range makes f^d negligible and we can keep only the term of fi^ucr when discussing the 
process of capture and escape of Cs atoms. 

It can be seen from formula (3) that fitter only relates to the diagonal elements of 
density matrix so we can use rate equations to calculate it. Taking into account the effect 
of spontaneous emission, stimulated emission and stimulated absorption, we have the 
following rate equations 1111 


■L'* = -j[s m , m+l (cr') 2 (py 


p { :i j + 
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and the normalization condition 


I 




(7) 
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After obtaining the steady solutions p { £ and p ( n f ) to the above equations we can 
calculate the scattering force exerted on a Cs atom under different physical conditions. 
Fig.(a) shows the influence of light intensity, Fig. 3(b) is the variation of force with different 
laser detuning;Fig.3(c) shows the forces as a function of the magnetic field gradient at a 
fixed point in space In these figures, the X-axis represents velocity in m/s and the Y-axis 
is the force in units of /ikr ( it equals to 4x 10' 21 N for the transition Cs 6S 1/2 ->6P 3/ 2 ). 




(a) For curves in the direction of the 
arrow, from lower to upper, 1=0.27, 
1.1 , 4.4 ,18 , 27 mW/cm 2 , at A - 
- 10.6 MHz, B=0T 


(b) For curves in the direction of the 
arrow, from right to left, A - - 10 6, 
-21.2 , -42.4 MHz, at I=4.4mW/cm 2 , 
B=0T 



(c) For curves in the direction of the 
arrow, from right to left , dB/dx=0.5 , 
1 , 1.5 mT/cm , at I = 4.4 mW/cin 2 , 
and A = -10.6 MHz. The distance 
to the center of the trap is 1cm 


Fig. 3 The force experienced by a Cs atom under different conditions 

Fig. 4 shows the position dependence of the scattering force on a Cs atom with 
different velocities The X-axis is the position and the Y-axis represents the force as in 
Fig.3 
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Fig.4 The scattering force of a Cs atom with 
different velocities . For curves in the 
direction of the arrow , from left to right 
V= -4.5, 0,4.5 m/s, at 1= l.lmW/cm 2 , 

A = -10.6MHz , dB/dx= 1 mT/cm . 


From Fig . 3(a) it can be seen if an atom is stopped at the center of the MOT and has a 
small velocity , we can obtain 

f ,ca,"r = - a V (8) 

From Fig 4 , we derive that if an atom has zero velocity and is near the trap center, 
then 


fscaner = ~kX (9) 

Therefore when both X and V are relatively small, we may rewrite f sclUer approximately 


as 


f'ca.ur = -kX - aV (10) 

where K and a are the string constant and damping coefficient respectively This is 
the typical form of force in a potential well with damping . 

B Calculation of the capture and escape of Cs atoms 

Before the calculation , we would like to explain several related concepts and introduce 
two important parameters . In practice the MOT is three dimensional and is formed at the 
intersection of six laser beams . The size of the MOT is determined by the radius of the 
laser beams . We will then also use the radius of laser beam as the spatial range of 1D- 
MOT The atom is considered to be trapped if it finally stops at the trap center The 
maximum velocity of atoms captured at the edge of the MOT is defined as the capture 
velocity V c . Mainly due to collisions with fast background atoms , trapped atoms at the 
center of the MOT may obtain enough initial velocity to be knocked out of the trap The 
minimum initial velocity is defined as the escape velocity V e 
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For the calculation of the parameters of motion the Cs atom is treated as classical 
particle . Its dynamical behavior obeys the Newton’s second law The force is obtained by 
the approach discussed above After integrating the equations of motion by the Runge- 
Kutta algorithm . we obtain curves showing variation of the velocity and position of a Cs 
atom in the process of capture and escape . Fig . 5(a) shows the variation of the position 
in the capture process . The motion of Cs atoms with different initial velocities is given 
where the atoms’ initial position is supposed to be at the edge of the trap Fig 5(b) shows 
the variation of the velocity in the same process . Finally Fig 6 presents the same curves 
in the escape process . Here the atoms’ initial position is assumed to be at the trap’s center 




(a) Variation of position (b) Variation of velocity 

Fig .5 The variation of position and velocity in the process of capture . For curves in the direction 

of the arrow , from upper to lower , the initial velocity is V, = 0 , 10 , 12.58 , 15 , 20 m/s , the 
radius of the laser beam is 2cm , I = 4.4 mW/cm 2 , A = - 10.6 MHz , dB/dx = lmT/cm 
From these curves , we obtain V c = 12.58 m/s . 



(a) Variation of position (b) Variation of velocity 

Fig .6 The variation of position and velocity in the process of escape . For curves in the direction 
of the arrow , from upper to lower , V, = 5 ,7.5,9.08 .10, 20 m/s . The parameters of laser 
are the same as those in Fig .5. From these curves we see that the escape velocity evaluates 
to be V e = 9.08 m/s . 

V c and V e are determined respectively in the following way : first , we set the velocity 
range to be [Om/s, 20m/s] , which , according to our calculation , cover the value of V c 
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and V, ; then , the velocity range is step by step decreased by the dichotomy , finally , we 
get the approximate values of V c and V e for a velocity range smaller than 0.01 m/s . The 
magnitude of V c and V e is related closely to the number of trapped atoms 


N 


s s 



[ 5 ] 


( 11 ) 


where r is the radius of the laser beam , n is the average velocity and a is the cross 
section for a trapped atom to be knocked out of the trap by a background atom.Obviously 
, o decreases with V e . From (11), we know that for a large number of trapped atoms it 
is beneficial to increase V c and V e by changing the trap parameters 

3 . DISCUSSION 


A . The properties of the force in MOT 

A Cs atom experiences three kinds of forces in MOT . These are the gravitational force 
the magnetostatic force and the optical force . The magnetostatic force is caused by the 
inhomogeneous magnetic field . When dB/dx = 2 mT/cm , a cs atom in the F g = 4 ,m g = 4 

sublevel obtains an acceleration of about g (gravitational acceleration ) The optical force 
is due to the interaction between the light and the atom . For A = - 10 6 MHz and I = 
4mW/cm , the acceleration caused by this force is about 103 g , which is much greater 
than that by the gravitational force and the magnetostatic force Therefore , we could 
neglect the gravitational and the magnetostatic force . 

As discussed above, both f,c lter and contribute to the optical force in the ID-MOT 
The capture and escape of atoms are mainly determined by the former due to its large 
interaction range of velocity and space Even though f^d has a smaller interaction range 
than that of fitter but it acts in this very small range much stronger on the atom [10] 
When atoms in MOT reach a steady state , they concentrate at the trap center with very 
small velocity amplitudes In this situation , the temperature and density of atomic cloud 
are determined by f^d 

B Contribution of Zeeman effect and Doppler effect 

In MOT the origin of f is due to the different relative detuning for <j* and o' 
photons This imbalance between the absorbed photons from opposite directions results in 
the net scattering force But an atom in the MQT can not tell whether the relative 
detuning is caused by the Zeeman effect or by the Doppler effect , because there is certain 
equivalence between these two effects on their contribution to the net force When 
Zeeman splitting between A m = 2 sublevels of a Cs atom excited state in compensated 

by Doppler shift , i.e. 

„ . UB 

gtMa ir x = kv (12) 
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then the relative detuning for cr* light is the same as that for a , Therefore the scattering 
rate of photons is equal in both directions , under which circumstance the atom 
experiences no force , 

/ = 0 (13) 

By substituting (12), (13) into (10) , we obtain 

i dB . .. 

K = g e P B (hk) 1 — • a (14) 

which indicates the equivalence between the Zeeman effect and the Doppler effect , i.e. K 
is proportional to a . It also shows that they both are dependent on the parameters of the 
laser field (detuning and intensity ) However they are different in that K is proportional 
to dB/dx whereas a is independent of dB/dx , which indicates that the scattering force 
originates in the inhomogeneous magnetic field . Under the condition of small V and X the 
calculation is simplified and we have obtained the approximate formula for a and K (see 
appendix for detail ) 


a ~ -0.9644 


hk 2 ss 0 r 


(1 +S 0 )(6 2 + r 2 74) 


K * -0.9644 


g<p B kss 0 r 


dB 


(1 + S 0 )(<5 2 + r 2 / 4) dx 


(15) 

(16) 


C . The relation of the force on the parameters of MOT 

We will concentrate our discussion on the effect of the laser intensity I , the frequency 
detuning A , the radius r of the laser beam and the gradient of the magnetic field dB/dx on 
the scattering force in MOT 

From Fig. 3(a) we know that for small intensities the force increases sharply with the 
laser intensity When 1 > I 0 S 2 T' 2 = 4.4m W / cm 2 ( i.e. S 0 = 1. The saturated intensity 
for the cycling transition F g = 4,m g = 4 — > F e = 5 ,m e =5 is I 0 = 1.1 mW / cm 2 ) , the 
intensity has only a small influence on the force S ’ m m 1 1 increases with the laser intensity 
but the population difference p l m s) — will become small due to the saturation effect 
Therefore the influence of laser intensity is diminished due to the compensation of these 
two effects with each other Fig 3(b) shows that the capture range of velocity increases 
with laser detuning Since the volume of the MOT* is determined by the radius r of the 
laser beam , increasing r may increase the volume efficiently ; Fig 3(c) shows the effect of 
the magnetic field gradient . At a point far from the trap center , atoms will be cooled to 
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non-zero velocity . The curve of the force is asymmetric about its zero point , which is 
caused by the difference of g-factor between the upper and the lower state . ( Upper state 
6P V2 ,F -5 ,g e = 0.4 , lower state 6S V2 ,F = 4 ,g g = 0.25 ) From Fig.3(c) we know that 

increasing the gradient dB/dx tends to decrease the damping force, whereas the string 
constant K increases with dB/dx(See [16]) This situation sets a certain limit on the value 
of dB/dx. 

D. Several measures to increase the number of trapped atoms 

We have known that increasing light intensity will not increase the number of trapped 
atoms if I>Io S 2 T ' 2 . Under the condition of I * Io 5 2 T ’ 2 , enlarging r (the radius of the 
laser beam) may increase the volume of MOT and increasing 6 ( laser detuning ) may 
increase the velocity capture range. After r and 8 being determined, we use the condition 
§m,m±i ^ 0 as the limit on the value of dB/dx. Among all the relative detunings , 5 ^.^(V = 
0,X = r) = 5 + ^' 1 rp B dB/dx is the biggest one. If it is negative, i.e., (note that 8 < 0) , 
the condition will be satisfied 

All in all, we think that the following methods may increase the number of trapped 
atoms: (1) increasing the detuning; (2) enlarging the radius of the laser beam as possible as 
to keep I « I 0 8 2 T 2 ; (3) dB/dx = -h 8(p B r)’ 1 when r and 8 are determined 

4. CONCLUSION 

The properties of MOT have been studied profoundly for years. On the basis of these 
studies, we used a simple method to calculate the force and motion of Cs atoms in 1D- 
MOT. The approximate formulas for a and K are given in the case of small velocity and 
volume and the characteristics of the force and the contribution of Zeeman effect and 
Doppler effect are also discussed The condition to increase V c and V e upon which we 
may build an efficient MOT are given 

However in this paper we don’t discuss the problems about the atomic density and 
equilibrium temperature in MOT, because the approximate method we used does not 
include f gr,d and atomic momentum diffusion. Thus it cannot be applied to calculate the 
atomic density and temperature. Furthermore the theory of density and temperature in 
MOT is quite complicated because it requires the consideration of various effects, which 
are not easy to evaluate. Anyway the related work is in proceeding. 

The authors are gratefully indebt to Prof. Wu Jike of Department of Mechanics, 
Peking University and Dr. Zhang Wenqing of Institute of Physics, Academia Science 
Sinica for their instruction and help with the computer program 


APPENDIX 

At first, we give the cited definition of the saturation factor 
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s n = 


q 2 / 2 
s 2 + r 2 / 4 


(a 1) 


and the relative detuning when X is small can be rewritten as follows 


S ^ 


(a.2) 


Finally the approximate population difference between the upper and the lower state is 


( g ) 


o ( ° 

P m + 1 


( 0 ) 


1 + 


(a 3) 


where /O ( m 0) is the population of sublevels of the ground state in the case of very weak light 
intensity under a + -a ' configuration and the calculated results are presented below 


< 0) = p ( 7 = 0.3 864 


>< 0) = = 0.0498 


p\ 0) = p = 0.03 3 6 , 

p { 0 0) = 0.0172 

After combining (3), (a.2) and (a 3), we have 
hkY . Q 2 / 2 
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= p w = 0 0188 


2(1 + S 0 ) (S - kV + r A 4 


Z pL o ;(c;;*’) : - 


Q 2 / 2 


(<5 + &K) +r /4 


)’i 


(a4) 


By substituting the values of p „ and corresponding C-G coefficients 


I p»> ( c;;-v = x P ': pc : : -')■ = 0.4 82 2 


(a 5) 


and using the approximate condition 5-kV = 6+kV « 5 for very small V as well as formula 
(a 1) , finally we obtain 
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scatter 


0 96 4 4 h_k 2 _SS ± V 
(1 + S 0 )(£ 2 + T 2 / 4) 
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a 


-0.9644 


hk 2 SS 0 r 

(1 + S 0 )(S 2 + r 2 /4) 


(a 7) 
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